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We study the internal dynamics of bosonic atoms in an optical lattice. Within the regime in 
which the atomic crystal is a Mott insulator with one atom per well, the atoms behave as localized 
spins which interact according to some spin Hamiltonian. The type of Hamiltonian (Heisenberg, 
Ising), and the sign of interactions may be tuned by changing the properties of the optical lattice, 
or applying external magnetic fields. When, on the other hand, the number of atoms per lattice site 
is unknown, we can still use the bosons to perform general quantum computation. 



I. INTRODUCTION 



In the last years there has been a great development of the storage and manipulation of cold atoms in optical 
lattices [lll3- The interest of these atomic systems is multiple. First, the optical lattices are accurately described by 

fTi the Bose-Hubbard Hamiltonian, and exhibit a quantum phase transition from a superfluid to a Mott-insulator ^. 

ri The jump from superfluid to insulator may be explored in the experiments in a time-dependent way jj, opening a 

wonderful playground for the study of time dependent phase transitions, establishment of coherence, and many other 
complex phenomena. Second, as it was shown in Refs. [a|, the bosons in the optical lattice may be used to simulate 

S ■ different spin Hamiltonians, with the advantage over magnetic materials that the parameters may be changed faster 
and more easily. Finally, these systems are also a good candidate for performing scalable quantum computations 



O 



'"^ , In this paper we focus on two applications of optical lattices: quantum simulation of spin s = 5 Hamiltonians 

^ ■ and universal quantum computing. Our proposals require tools which already belong to current experiments: (i) an 

O _ optical lattice with one or more atoms per site, (ii) atoms with two accessible hyperfine levels, (iii) Raman lasers that 

connect these levels, and (iv) a magnetic field or an electric field with the appropriate spatial dependence. With these 

tools, and assuming that the lattice is a Mott-insulator with one atom per lattice site, we generalize the techniques 

developed in Ref. 5i] to consider also these external elements (Rabi oscillations and magnetic and electric fields), to 

^ ' engineer a large class of spin Hamiltonians, 
00 ■ 

^ : H = Y, [Xl^ata] + \t, {afa^ + af aj)] + ^ ks.. (1) 

S ■ ^'^'^ 

^P \ The flexibility of the system is such, that one may simulate an Ising or Heisenberg Hamiltonian, on different geometries 
^^ ■ dictated by the underlying optical lattice, and even with some randomness in the coefficients \z,± or in the effective 
magnetic field, h. 

C^ ■ The second focus of this work is quantum computation. The effective spin Hamiltonians mentioned above could 

be used to develop a universal quantum computer, either directly, or by constructing the so called cluster states |3|. 
However, since it is not easy to produce a bosonic crystal with integer filling factor — or at least not with the fidelity 
' O ' required for scalability — , we have developed other methods for quantum computation with any number of atoms per 

C , lattice site. The proposal that we developed in 9] and summarize here combines ideas from quantum engineering (for 
the definition of qubit and to obtain an effective interaction between them), the technique of adiabatic passage [lO| 
(to perform controlled gates even when the parameters of our system can not be accurately determined) , and ideas 
from quantum control [ll| (to cancel dynamical phases in the adiabatic passage). 

The paper is structured in two parts. In the first part we concentrate on the dynamics of optical lattices in the 

5_^ • Mott insulator regime. The lattice will be loaded with one atom per site, and each atom has two accessible internal 

d I state, so the atoms may be formally identified with s = \ spins. We will write down the most general Hamiltonian 
that takes into account the hopping and interaction of atoms, as well as the influence of external magnetic and electric 
flelds on the atoms. By adjusting the optical lattice and the external fields, we will be able to develop an effective 
Hamiltonian in the spins representation which covers all possibilities shown in Eq. (^. We determine the accuracy 
of this description analytically and numerically, and brieffy study the application of these techniques to generate 
entanglement between the atoms. In the second part of the paper we will demonstrate that, if there are imperfections 
in the loading of the lattice, the Hamiltonian of the system is not known with enough accuracy to do quantum 
computing in a "traditional" way. Next we develop a technique to circumvent our ignorance about the Hamiltonian, 
using adiabatic passage with the different parameters of our Hamiltonian to produce a universal set of gates. Finally 
we estimate the errors of our proposal, studying both the inffuence of the speed of the adiabatic process, and of 
imperfections in the setup. 



II. SIMULATION OF SPIN HAMILTONIANS 
A. The Bose-Hubbard model 

In an optical lattice, pairs of laser beams create a stationary wave which the atoms see as a periodic potential. If 
the energies involved in the dynamics are so small that the second Bloch band never gets populated, we may use the 
Bose-Hubbard Hamiltonian to describe the atomic ensemble Q • We will assume that the lattice is populated with a 
single atomic species, and that only two of the hyperfine levels of these atoms may be excited. Then the Hamiltonian 
reads 

4o? = ~Ja{a\a,+h.c.) + Mh\b,+h.c.), (2b) 

-ffmt = \Uaaa\a\a,a, + \Ut,bb\b\bA + Uaba\h%a, (2c) 

Hl^ig = ^k{a\a,-b\h{), (2d) 

H^i = lMa^ + b\h), (2e) 

^E = ^{Ahe^^ + b\a,e-^^). (2f) 

The operators a and b are bosonic destruction operators for atoms in two degenerate hyperfine states; the indices 
i and j run over the lattice sites, and (i, j) denotes a pair of neighboring sites. The constants J and U depend on 
the depth of the optical lattice, and they measure the amplitude of probability of atoms hopping to neighboring sites 
(Hhop) and their effective on-site interaction (Hint), respectively. The Hamiltonians Hmag and Hei account for all 
possible energy shifts on the internal states of the atoms. They may be generated by means of magnetic fields (a and 
b represent two different hyperfine states which suffer different Zeeman shifts), or with highly detuned laser beams 
which induce a Stark shift. Finally, Hias models transitions between the two internal states of the atoms. The Rabi 
frequency fi and the phase (j) are related to the intensity and the phase of the laser which induces these transitions. 

Throughout the paper we will assume that the system is in the Mott-insulator regime, U ^ J, in which the hopping 
represents a small perturbation with respect to all other term of the Hamiltonian. For the external fields, e^ and 7^, 
we will require that they have a simple spatial dependence, i. e. they may increase or decrease linearly along a given 
spatial direction, as in e/c = cq -I- (5 x k. No particular boundary conditions are imposed, and all results may be trivially 
generalized to any geometry of the optical lattice. 

It is worth remarking here that, nowadays in most experiments there exists a residual harmonic potential which is 
used to further confine the atoms within the optical lattice. If the gradient of this potential is extremely small compared 
to J and U, it may be regarded as a spatially dependent contribution to efe. In general, however, the harmonic 
confinement influences greatly the ground state of the Hamiltonian ^ creating coexisting regions of superfluid and 
insulator phases [3| which are useless for our purpose. It would be possible to get rid of this potential if we use 
additional optical elements to create a barrier that prevents the atoms from escaping through the borders of a lattice. 

B. Effective spin Hamiltonians 

It is well known that for a Mott insulator made of one bosonic species the Bose-Hubbard model is equivalent to 
the XY model, where spin states are identified with holes and particles in the lattice jj. More precisely, if the filling 
factor |14j is in between n and n -I- 1 and the on-site interaction is strong, then the occupation number of each site will 
be either |n) or |n -I- 1). We identify these states with spin states | + ^) and | — 5), to obtain the effective Hamiltonian 

He = -A± E«^i + ^^J) + E(^^ + ^^)<- (3) 

The magnetic interaction, Aj^ = J{n -1-1), is provided directly by the hopping terms, Hhop, and it can be rather 
intense. However, from the point of view of quantum simulation this approach has an important restriction: since the 
number of particles is conserved, the total spin of the system is constant, and we cannot introduce terms which are 
proportional to a^ or a^ in the Hamiltonian. 

We can simulate a larger family of spin Hamiltonians if we populate the lattice with a single atomic species 
with two degenerate internal states. We will assume that the lattice is loaded with exactly one atom per site, and 



identify the two possible states of each site, a|, \vac) and 6j. \vac), with the two polarizations of the spin, ji^)^- 
The states with single occupation are separated by an energy gap of order 0(U) from any other configuration of the 
lattice. Since we are deep in the Mott insulator regime, J <^ U, states with double or higher occupation are highly 
improbable, and we may treat the hopping term. Hi — H^op, as a perturbation with respect to other contributions, 
Hq = Hint + Hmag + H^i + Hias ■ Usiug sccoud Order perturbation theory ,12] we will write an effective Hamiltonian 
within the spin space, which looks as follows 

m,\j) = mo + Hi\j) + i ^ (z|i/iio \^^ + Y^] ^^i^^i-^'^ ■ (^) 

t: L * s J s - 

While the indices i,j run over the Hilbert space of the spins, ^ represents any configuration with an excess or defect 
of particles in any well. The numbers Ei, Ej and E^ are the unperturbed energies of these states: Ei = {i\Ho\i) , etc. 
We will analyze separately how the different terms in Eq. |(2J) influence the effective Hamiltonian |0J. First of all, 
if there are no external fields (e^ = 7^ = 51 = 0), iJg is a Heisenberg Hamiltonian (Q with constants given by Q 

1. _ 'J a + -Jb _ "Ja _ "^b \ _ _fJaJb_ 

2Uab Uaa Ubb Uab 

If we switch on a linearly growing electric field H^i (|2ell with 7^ 
same, but the constants change: 

, _ {Ja+ JlWab JlUaa J^Ubb ^ _ 
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■^{ui,-5^) m,-5^) {ui-5^y 

Effectively, the application of the electric field H^i is equivalent to a change of the interaction constants, of the form 

C/2 -,52 
Uuv~-^j , u,v€{aM- (7) 

This effect may be used to intensify the interactions, making the spins evolve faster. But it also may be used to 
introduce some randomness in the system. For instance, if 7^ does not grow linearly, but fluctuates from site to site, 
the constants, {A^, Aj^, K\ will also fluctuate from site to site, with an expression given by Eq. © with 5 -^ 7i+i — "it- 
Another interesting effect is provided by a linearly growing magnetic field Ij2d|l . tk = ^q + 5 x k. This contribution 
to the Hamiltonian breaks the degeneracy between the states | + 5 , — 5 ) and | — 5 , + 5 ) ■ If the gradient of the magnetic 
field is weak compared to the interaction {J ^ 5 -^ U), the system is still described by a Heisenberg interaction ^, 
with constants given by Eq. ©, except for the magnetic field 



JbUbb JgUg 

Ul-S^ ~ Ul.-6^- 



hz.i - ei + y.^2 _ p IJ2 X2 ■ (^) 



If the gradient of the magnetic field is comparable to the interaction {S ~ U/20, for instance), the splitting between 
l+i, — i) and I — i, +i) becomes so large, that they may no longer be connected by the Hamiltonian. Applying the 
rotating wave approximation one finds an Ising Hamiltonian, 

He^Y. [>^-^t^Ul + K^(^t] ^ (9) 

i 

where X^ is still given by Eq. ©, and hi by Eq. |(SJ). Another way to achieve an Ising Hamiltonian which was shown 
in Ref. |5j consists in tuning the properties of the optical lattice so that Jb -^ 0. This way we set A_l = in Eq. ^, 
but there is a residual magnetic field, h^ = —J^/Uaa- 

All Hamiltonians which we have shown up to now commute with the operator Sz — Yl '^t ^ ^a ~ Nb- In order to 
change this component of the spin, we have to allow transitions between the internal states a and b. Such processes 
are modeled by the Hamiltonian Hias (Efl), and with little work it can be shown that these terms translate into an 
effective magnetic field along the X and Y directions of the spin, J^ii^xO'f + hyof), with 

/ix = -^ cos(0), ^y = -^ sin(0). (10) 

This last term completes all what is required to simulate the family of Hamiltonians presented in the introduction 
(|T|l. We may now particularize the previous results for current experiments with Rubidium. If we do not use state- 
dependent optical lattices (i. e. all atoms see the same potential), we can approximately take Uaa = Ubb = Uab = U , 



and Ja ^ Jb ~ J- Then, the effective Hamiltonian with electric field Eq. QJ becomes 






and the Ising Hamiltonian with the linearly growing magnetic field is written as 



He = Y. 



(■i^z - 



C/2-J2 



(12) 



In this case, for simplicity, we have omitted corrections which are proportional to j'^S/iU'^ — (5^), but which only take 
place at borders of the lattice. 

Of all Hamiltonians which we have shown, the Ising interaction has the greatest interest for quantum computation, 
because of its simple form and because it produces a universal gate known as the controUed-phase gate jis]. However, 
for many purposes it would be interesting to get rid of the effective magnetic field which appears both in Eqs. (jS)) 
and (|12|l . and which introduces an uncontrolled dephasing. We know two ways to produce an Ising Hamiltonian 
without magnetic field. The first one would be to use Feschbach resonances to increase the interaction between atoms 
of type a and h, while keeping equal tunneling rates. In other words Uaa,Ubb ^ Uab and Ja = Jb- If this is the 
case, then both Aj_ and h become zero in Eq. ^ and we get the desired model. However, Feschbach resonances 
require intense magnetic fields which may be difficult to stabilize for a time scale of order Uaa.bb/ Jab- '^^'^ other 
method to get rid of the effective magnetic field is to set up the conditions which lead to Eq. 112() and then apply a 
spin-echo 7r-pulse at times T/2 and T, where T is the total duration of the experiment. This way the fc-th lattice site 
acquires a phase exp(— jft,fcCT^T/2) during the first half of the experiment, which is canceled with the unitary operation 
a% exp(— i/ifecr|T/2)crf of the second half. 

C. Error bounds 

In deducing the effective models (Q, 111|) . and (|12|l we have performed several approximations. For the Ising 
model, the first source of error relates the contributions of ct^ct^"'"^ and ayay'^^ . These operators induce a swap 
between qubits, so that the state | + ^, — ^) is connected to the | — 5, +^) and viceversa. Using the interaction picture 
and the Born approximation, we can estimate the probability of a swap happening in a pair of neighboring wells as 



pfc 



t 

2<5 



dre''' (+ - \V\ - +) 



- (52([/2-,52)2- ^ ^' 



Another possible source of errors which affects all models is the accumulation of particles on one well. Put in other 
words, the Hamiltonian ^ allows processes in which one particle jumps from one well to a neighboring one and 
remains there. To first order in perturbation theory, an upper bound for this probability is 






dre'^^-^^ {i\V\0 







' 2J2 



This probability is larger than the one for qubit swaping (|13|l . which means that in the worst case, the error when 
considering M wells is at most of order E = MO [j'^/{U — (5)^). 

The last and most dangerous source of error is the possibility of having in the initial state a cell with more than one 
atom. If we are simulating the Heisenberg model, these extra atoms may hop to neighboring sites. The scattering of 
these atoms with those of opposite polarization leads to very fast changes on the lattice, and our picture of localized 
spins breaks down. Therefore it is extremely important in these experiments to begin with a well prepared Mott 
phase with filling factor n = \. 

D. Fidelity 

We have applied our model Hamiltonian (|12|l to study how entanglement is created between two, three and more 
atoms in an optical lattice. The idea is not only to perform numerical experiments which show us how much entan- 
glement can be produced in a realistic experiment, but also to check the accuracy of our effective Hamiltonian ([T|) 
when describing the process. 




Co 



1 


-r 


^^\ 








(b) 






t 


\\ 


-^ 


0.8 


- / 


\ 


\ 


1 




/; 


0.6 


- / 






1 


\\ 


1 

1 - 
\ /' 


0.4 


-/ 




\i 






\ /' 

V /' " 

V /' 


0.2 


/ 




V 






V // 
W; - 



X 10 



500 1000 

TU 



1500 



FIG. 1: Evolution of entanglement as a function of time for an initial state given by Eq. H16|l . evolving under the conditions 
which lead to the effective Hamiltonian of Eq. Q . (a) We plot the entropy of the reduced density matrix of one well (solid 
line) out of two, when jumping from J — 0, 5 = ek+i — e^ ~ to J = O.OIL^, 5 = 0.05C/. In dash-dot line we also show the 
fidelity of the operation with respect to the effective Ising Hamiltonian H12|l . (b) We perform a similar plot for a system with 
four (dashed line) and five (solid line) lattice sites, for J = 0.0411 and S = 0.3C/. In dash-dot line we plot the fidelity of the 
effective Hamiltonian (1121 when describing this evolution. 



Our numerical experiments typically proceed as follows. First we prepare an initial state which has one atom per 
lattice site, and which we can identify as a spin state. Then we choose the parameters of the optical lattice (0) so that 
the theory developed in Sec. Ill Bl applies, and simulate the evolution of the whole system considering also the states 
with double and higher occupation numbers. Finally we measure the fidelity of the effective Hamiltonian, H^, as 



T 



\{iPiO)\e'"'e' 



iH^t 



1^(0)) 



(15) 



where H is the full model depicted in Eq. ^. The deviation of this quantity from 1 measures the error committed 
because of trying to describe the evolution with our simplified model Q) . 

We have performed two sets of simulations in which we check the fidelity of the Ising and of the Heisenberg 
Hamiltonians, respectively. In the first row of numerical experiments the starting point is a product of each atom 
being either in state a or 6 



1^(0)) = 



M 



U'=l 



M 



k=l 



M\+\) + \-\)) 



■M). 



(16) 



From a physical point of view, this configuration may be prepared by loading the lattice with one atom per site, 
switching the tunneling off, J ~ 0, pumping all atoms into internal state a, and finally applying uniformly a 7r/2 laser 
pulse over the whole lattice. Beginning with the previous state we customize the optical lattice so that the effective 
interaction corresponds to the Ising model (|12|l . For the results in Fig. ^a) we have lowered the optical lattice 
until J = O.Olt/ and set up a magnetic field with 5 = 0.05C/. As shown in Ref. [g, the ensemble of atoms should 
evolve periodically from a product state to an entangled state with maximal connectivity and large entanglement 
persistence, also known as clusters. To measure the entanglement of one qubit with respect to the rest, we compute 
the von Neumann entropy S{pk) — Yl "■^k log2(Afc), where pk = J2k ^k IV'fe) (V'fcl is the reduced density matrix of one 
well. In Figure ^a) we show how this value evolves for a system of two wells and two atoms. In order to make the 
process faster, we have made numerical experiments with larger hopping strengths, J = O.OAU, and stronger field 
gradients, S — 0.3C/, for lattices with four and five wells. The results are plotted in Figure ^b). The maximum 
entropy is reached in all wells simultaneously for state which is close to a cluster. But now we notice that the process 
is not perfect. First, the due to the greater intensity of the hopping, the system gets a larger contribution of states 
with occupation n ^ 1. And second, the presence of these states alters the dynamics, so that for long times the 
fidelity decreases. 

Another method to create entanglement in the optical lattice would be not to use the Ising interaction, but 
to find an initial state which, under the influence of the effective Hamiltonian (|ll|l . leads to a larger amount 
of entanglement. The initial state we propose is a set of alternating spins \ip) — |01...0101). For two wells, 
after a time T = 7r[//8J^, we get the Bell state \4'~)- For three wells, the states we get are of the form 
\ip{t)) ex [(1 -I- 2e*6*) |010) + (1 - e^**)(|100) + |001))] /3. For more than three wells, we have calculated the states 




FIG. 2: Entropy for first, second and tfiird well as a function of time (solid lines), for a system with five particles, five wells 
and open boundary conditions, with J = 0.04 U and no external fields. We also plot the fidelity (dashed line) with respect to 
the Heisenberg Hamiltonian lllll . 

numerically. As we show in Fig. |51 for open boundary conditions the entropy is not uniformly distributed and a 
maximum cannot be reached simultaneously on all wells. Nevertheless, for a given fidelity, the entanglement of one 
well with respect to the others grows faster with the Heisenberg model than with the Ising model. 

III. QUANTUM SIMULATION AND COMPUTING WITH MORE ATOMS PER WELL 

In the previous sections we showed how to simulate spin Hamiltonians using cold atoms in an optical lattice and 
some external fields. The whole procedure relied on the assumption that one can load the lattices with precisely one 
atom per site. However, this is not the situation of current experiments, where the number of atoms per cell is small 
(1 to 3), but otherwise ignored. Motivated by this limitation, we have developed a scheme for quantum computation 
which works independently of how atoms are distributed over the lattice. As a possible application, this general 
quantum computer might be used to simulate spin Hamiltonians. 

The method which we summarize in this section has three key ingredients. First, we have to choose a set of states 
which are suitable for quantum computation, called qubits. Several conditions have to be imposed on the system 
to ensure that the different occupation numbers do not influence the computation. Second, we have to design an 
effective Hamiltonian which allows both a predictable interaction between qubits and the possibility to perform local 
operations on the qubits by means of external fields. At this point one realizes that if we do not know how the 
particles are distributed over the lattice, then we also lack a lot of information about the Hamiltonians. This is where 
the third ingredient comes in: we design a procedure to perform controlled unitary operations on the lattice, even 
when most of the parameters of our system are unknown. 



Basic tools 



To do quantum computing with m qubits, we need to restrict the possible configurations of our system to a 2™- 
dimensional Hilbert space, which is our computation space. Formally, this space is isomorphic to the space of ^-spins 
developed in Sect. IIIBI but we cannot use those definitions because they required to have only one atom per lattice 
site. Instead, our implementation of a qubit will rely on the number of atoms which have been excited from one 
internal state, a, to the other, b, 



|0) (X (a^)" I vac) , |1) ex 6t(at)n-i \yac) 



(17) 



In order to ensure the stability of the computation space, we have to impose certain conditions on our system. First, 
in the absence of external fields, the states |0) and |1) should have the same energy. Otherwise, different configurations 
of the lattice will acquire different phases which spoil the computation. To avoid these dephasings, we just need that 
Uaa = Uab- However, for convenience, we will first assume that Uaa and Uab are zero, and only at the end study what 
happens when this is not true. 

The second requirement is that our lattice remains in the computation space at all times. That is, there can 
only be at most one atom of type b per site, and the number of atoms per lattice site must remain constant. To 



avoid exciting more than one atom, there must exist a quantum blockade mechanism which separates the states 
\m) = b' a^ \vac) , m — 2,3 . . ., from our computation space and depopulates them. This blockade is achieved 

for Ubb 3> Uab, Uaa- And finally, to avoid problems with particles moving from site to site, we will impose a gradient 
of energy along the lattice as in Eq. if^ with ek — eo + S x k. The gradient 6 must be large compared to Jb and 
the possible residual values of Ja, Uaa and Uab, but it also has to be smaller than Ubb to prevent the actual motion of 
atoms. In other words, JbM, Uaa, Uab ^ \Ubb — (5| <C (5, Ubb- Under the previous conditions, an adiabatic elimination of 
the hopping term converts Eq. (0) into an Ising model 

Apart from an effective interaction between qubits, we also need means to modify the state of a certain lattice site. 
To do this we will operate on the qubits using light or magnetic fields. The interaction between the atoms and these 
external fields is governed by the Hamiltonians (|2d|) and l(^ 

H^j (4«fe - bib,) + ^ [albke^* + bla.e-^*) . (19) 

When we perforin an adiabatic elimination of the states outside our computation space, Eq. H19(l turns into 

H^ = ^a, + V^^ {a+e^^ + a_e-**) . (20) 

This effective Hamiltonian is problematic, because it depends on the occupation number, which is unknown. It is 
not be possible for instance to use a 7r/2 pulse to prepare the state (|0) + |l))/\/2 because, since do not know n, we 
also cannot predict the time during which the laser should operate. Furthermore, one might also argue that, due to 
the characteristics of the optical lattice we cannot control the parameters Jb, Ub and S with enough precision to do 
traditional quantum computing with it. 

All previous problems may be formulated in a more general way: we want to perforin accurately certain unitary 
operations on our set of spins or qubits using the Hamiltonians Hi ()20|) and 

^2 -y 111) (111, (21) 

even when we do not know some of the constants in them. To be precise, in what follows we will assume that A, 51, 
and A are unknown, but they can be set to zero and reach a positive value (A™, flm, and Am). The only parameter 
which will be precisely controlled is the phase of the laser, 0. 

There exists a solution to this problem which is based on adiabatic passage. According to the adiabatic theorem, if 
we change the parameters of a Hamiltonian slow enough, and the Hamiltonian has no degenerate eigenstates, we will be 
able to perform the unitary transformation U{T) = ^^ e^'^'^"^^"'' \^ a{T)) {^ a{0)\ , where |$Q,(i)) are the instantaneous 
eigenstates of our system. If the process is designed carefully, and the geometrical and dynamical phases are properly 
canceled {(j)a = 4'a = 0)i the resulting transformation does not depend on the precise values of the parameteis which 
governed the evolution (A, fl, etc), but on the path which the system followed in the space of possible Hamiltonians. 
The design of these paths will be discussed in the following section. 

B. Design of the quantum gates 

It is a well known result [ltj| , that a quantum computer can be built upon a small set of unitary operations. By 
combining these "gates" we may approximate any other transformation as precisely as required. In this subsection 
we will explain how to use the atoms in the optical lattice to perforin a universal set of gates made of a phase gate, 
U = e'^'^-/2, a Hadamard gate and a CNOT gate. 

First we care about the phase gate. In (|20|l we set A = for all times and change the remaining parameters {il, ip) 
as depictured in figure O^a): 

(0,0) ^ {nrn,o)^^ inm,0/2)^'^\nrn,e/2 + 7:) 

^^^ {n,n,e + 7r)^-^{o,e + TT) (22) 

All steps are followed adiabatically and require a total time T, except for step (iii) whose double arrow indicates a 
sudden change of parameters. Note that 17(0) = n{2T) = 0, ^(t) = n{2T-t) and ip{t) = 7r + 9-ip{2T-t), which does 
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FIG. 3: Schema of how the parameters of Hamihonian 12011 have to be changed in order to perform a phase gate (a), and 
Hadamard gate (b). 

not require the knowledge of the function / but implies a precise control of the phase. A simple analysis shows that 
(i-v) achieve the desired transformation |0) — > e'^^^ |0), |1) -^ e~**/^ |1), with a total cancelation of all geometrical 
and dynamical phases. 

The Hadamard gate can be realized in a similar way. In the space [A, n^ — ^ cos{ip)], the protocol is 



(0, n. 



i^) 



\^7n-) ^ ^71 



^^ (A 0) ^ (A n 



^^'(o,a„)^(o,-a„y"^A,-a: 



{vii) 



(A,0), 



(23) 



as shown in figure Ofb). In order to avoid the dynamical phases, we have to make sure that steps (i-v) are run 
in half the time as (vi-vii). More precisely, if i < T, we must ensure that A{t) = A(T — i), flx{t) = flx{T — t), 
A(T + i) = A(t/2) and fl^iT + t) = nx[t/2). With this requisite we get ^(|0) + |1)) ^ |0), ^(|0) - |1)) ^ - |1). 
Again, the whole procedure does not require us to know 51 or A, but rather to control the evolution of the experimental 
parameters which determine them. 

For the last universal gate we will take a pair of interacting qubits, and shine a laser on one of them. The total 
effective Hamiltonian may be written as 



H^^\ll) (111 +1® ^{^+e'^ + cT-e"'-^) 



By changing the parameters [A, 0.^ = ^ cos{(f)] in the following way. 



(A™,0)^(A 



„,51„) ^^ (0,r!„) ^ {0,-^m) ^^ (A. 



-51„)^(A„,0). 



we obtain the transformation 



t/i = |0)(0|®I+e*«|l)(l| 



I Idy 



(24) 



(25) 



(26) 



where ^ = /q S{t)dt is an unknown dynamical phase. To get rid of this phase we have to apply two more unitaries. 
First we need to perform a NOT on the first qubit U2 = (|0) (1| + |1) (0|) (g)I. Finally, if A(^^(i) denotes the evolution 
of A in equation (|2SJ), the system must follow a path such that A^^^t) = A^^^t), Q^^^t) = 0. If the timing is 
correct, we achieve U3 — (|0) (0| + e*^ |1) (1|) (g)X Everything combined gives us the CNOT up to a global phase 
Ucnot ~ |0) (0| (K)1I+ |1) (1| <E) icTy = e"*^ [72^3^/2^1, which does not affect the computation. 



Errors 



In our design for a quantum computer, there are many implicit approximations, which in practice become sources 
of error. Basically, we have neglected processes which lead to (i) excitation of more than one atom into state |6), 
(ii) change of occupation numbers due to hopping of atoms, and (iii) hopping and interaction of atoms in state \a). 

The first two phenomena are suppressed if {Q/Ubb)^ ^ 1 and {Jj^ /Ubbf ^ 1- We may analyze the remaining errors 
perturbatively, and study how they modify the effective Hamiltonians H2()|l and (|21|l . First, the the virtual excitation of 
two atoms increments the parameter A by an unknown amount, Ag// ~ A + 2VL^nk/{A + Uab — Ubb)- But \iUab ^ Ubb 
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FIG. 4: Log-log plot of the gate error, E = 1 — T, for the (a) phase, (b) Hadamard and (c) CNOT gates. The parameters for 
the simulations are Uaa = Uab = Jm, Ja = Jb,g = Ubt + Uab/'i, and T — lOO/fim. For the local gates we choose Am ~ 6Slm = 1, 
and for the nonlocal gate Qm = Jm/&- For each simulation we choose different population imbalance (|n — »n| = 0, 1, 2 for solid, 
dashed and dotted lines), and change the interaction constant Ubb- 

and VPrikT /Ubb ^ 1, this shift may be neglected. And second, in the two-qubit Hamiltonian (|21|l there appear 
additional contributions due to virtual hopping of all types of atoms, which are of order niax(J(,, JaY /S"^ ^ J"^ /Ubb- 
For J^T/Ubb ^ 1 these energy shifts may also be neglected. 

To check the validity of our approximations, we have simulated the evolution of two atomic ensembles with an 
effective Hamiltonian which results of applying second order perturbation theory to Eg. Q, and which takes into 
account all important processes. As a figure of merit we have chosen the gate fidelity [l^ 



J^=2-^\TT{Ul^^lUreal}\' 



(27) 



where n is the number of qubits involved in the gate, Uideai is the gate that we wish to produce and Ureal is the actual 
operation performed. The results are plotted in Fig. 0| In these pictures we show the error of the gates for simulations 
in which all parameters are fixed, except for Ubb and the occupation numbers of the wells. The first conclusion is that 
the stronger the interaction between atoms in state \b), the smaller the energy shifts. This was already evident from 
our analytical estimates, because all errors are proportional to 1/Ubb- Typically, a ratio Ubb = lO^'C/afc is required to 
make J- = I — 10^^, but reasonable fidelities may be achieved for more realistic values. Finally, the larger the number 
of atoms per well, the poorer the fidelity of the local gates [Figure 01 a-b)] . 



IV. CONCLUSIONS 



We have shown that a Mott insulator of bosonic atoms with filling factor one and two accessible internal states 
behaves as a lattice of localized spins. This result generalizes that of 5], by considering not only atoms in an optical 
lattice, but also the infiuence of external fields, such as lasers, magnetic fields, Rabi oscillations, etc. As we showed, 
tuning the parameters of the potential which confines the atoms and using additional elements, such as a magnetic 
field or lasers, the spins may be forced to simulate a great variety of Hamiltonians, including the Heisenberg and Ising 
models, with either attractive {Uab > 0) or repulsive {Uab < 0) interactions. When, as in current experiments, the 
number of atoms per lattice site is unknown — for instance there is some superfluid component, the filling factor is not 
an integer or there is some external potential — , we may still draw the atoms into a regime in which they behave as 
localized spins, but we will not know the intensity of the interactions of the system with external fields. Nevertheless, 
even ignoring the parameters which govern their dynamics, we have shown that it is possible to perform a universal 
set of gates and do quantum computation. 
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